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ABSTRACT 

Nucleolin is a multi-functional nucleolar protein that 
is required for ribosomal RNA gene (rRNA) tran- 
scription in vivo, but the mechanism by which 
nucleolin modulates RNA polymerase I (RNAPI) tran- 
scription is not well understood. Nucleolin depletion 
results in an increase in the heterochromatin mark 
H3K9me2 and a decrease in H4K12Ac and H3K4me3 
euchromatin histone marks in rRNA genes. ChlP- 
seq experiments identified an enrichment of 
nucleolin in the ribosomal DNA (rDNA) coding and 
promoter region. Nucleolin is preferentially associ- 
ated with unmethylated rRNA genes and its deple- 
tion leads to the accumulation of RNAPI at the 
beginning of the transcription unit and a decrease 
in UBF along the coding and promoter regions. 
Nucleolin is able to affect the binding of transcrip- 
tion termination factor-1 on the promoter-proximal 
terminator TO, thus inhibiting the recruitment of TIP5 
and HDAC1 and the establishment of a repressive 
heterochromatin state. These results reveal the 
importance of nucleolin for the maintenance of the 
euchromatin state and transcription elongation of 
rDNA. 



INTRODUCTION 

Nucleoli play essential roles in ribosome biogenesis and 
particularly in the transcription of ribosomal RNA 
(rRNA) genes. In a typical human cell, there are ~400 
copies of ribosomal DNA (rDNA) arranged in tandem 



repeats in nucleolar organizer regions (NORs). Each 
rDNA encodes a precursor transcript (45 S pre-rRNA) 
that can be processed to generate 18S, 5.8S and 28S 
rRNA. The sequences that encode 45S pre-rRNA are 
separated by long intergenic spacers (IGSs). 

rRNA genes are transcribed by RNA polymerase I 
(RNAPI) and require several associated factors for tran- 
scription specificity and efficiency. Recent works have 
highlighted the importance of the epigenetic control of 
rRNA transcription (1-5). In eukaryotic cells, only 
about half of the rRNA genes are transcriptionally 
active and can be distinguished by different epigenetic 
marks such as rDNA methylation, histone post-transla- 
tional modifications and different chromatin-associated 
proteins (6). Actively transcribed rRNA genes possess eu- 
chromatin marks characterized by H4K12 acetylation, 
tri-methylation of H3K4 and hypo-methylated CpG 
dimers, while silent rRNA genes have heterochromatin 
features such as di-methylated H3K9, tri-methylated 
H4K20, tri-methylated H3K27 and hyper-methylated 
CpG dimers (5,6). Apart from the post-translational 
modifications of histones, specific incorporation of 
histone variants such as H2A.Z and macroH2A (7,8) 
may also contribute to the epigenetic control of rRNA 
transcription. Araya et al. (2010) recently observed the 
differential incorporation of macroH2Al and two 
subtypes into the ribosomal cistron in fish could regulate 
gene expression during the acclimatization process (7). 

A major mechanism for silencing active rRNA genes 
seems to be the recruitment of the nucleolar remodelling 
complex (NoRC) to the promoter by the transcription 
termination factor-1 (TTF-1), which binds to the 
promoter-proximal terminator TO (1,9,10). NoRC, a 
member of ATP-dependent chromatin remodelling 
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machines, comprises the ATPase SNF2H and TIP5. Once 
NoRC is recruited to the rDNA promoter by TTF-1, it 
mediates heterochromatin formation and silencing of 
rRNA transcription by recruiting histone-modifying 
enzymes histone deacetylases (HDAC), histone methyl- 
transferases (HMT), DNA methyltransferase and by 
shifting the promoter-bound nucleosome into a silent 
position (1,2,10-12). Recent reports also showed that 
pRNA, a non-coding promoter-associated RNA, could 
form a triplex structure with TO, leading to displacement 
of TTF-1 from TO. The triplex could then recruit 
DNMT3b to the rDNA promoter, thus methylating 
CpG-133 and contributing to the repression of transcrip- 
tion (5,13). NoRC-dependent rDNA silencing and hetero- 
chromatin formation has been studied in detail. However, 
little is known about the mechanisms that counteract het- 
erochromatin formation and promote the establishment 
and maintenance of the euchromatic state of active 
rDNA repeats. 

Since it was first described as one of the main nucleolar 
proteins, nucleolin has been shown to be implicated in 
many steps of ribosome biogenesis including the synthesis 
of rRNA by RNAPI (14-20). Multiple functional domains 
allow the interaction of nucleolin with numerous proteins 
and nucleic acid sequences (RNAs and DNAs). Earlier 
experiments suggested a link between the proteolysis of 
nucleolin and RNAPI transcription elongation (21) and 
that only RNAPI (not Pol II or Pol III) transcription 
could be regulated by nucleolin, independent of the 
sequence of the transcribed RNA (15). Nucleolin deple- 
tion in different cell lines using small interfering RNA 
(siRNA) (14,19) or by conditional knockout in DT40 
cells (20) results in the reduction of the accumulation of 
pre-rRNA. Metabolic labelling and analysis of the matur- 
ation or pre-rRNA produced in absence of nucleolin 
strongly suggest that nucleolin is required for efficient 
transcription of rRNA genes (20). Although the mechan- 
ism of nucleolin action on the synthesis of pre-rRNA 
remains unclear, several experiments indicate that this 
regulation may be achieved through chromatin. 
Nucleolin binds tightly to chromatin (22,23) and is able 
to modulate chromatin structure by interaction with 
histone HI (24,25) or to stimulate the remodelling 
activities of the ATP-dependent remodelling complexes 
SWI/SNF and ACF on canonical or macroH2Al contain- 
ing nucleosomes (26). In vitro, nucleolin possesses a 
histone chaperone activity. It binds directly to H2A- 
H2B dimers and facilitates their assembly into nucleo- 
somes on naked DNA (26). FRAP experiments on 
eGFP-tagged histones (H2B, H4 and macroH2A) in 
living cells depleted in nucleolin also confirmed a role 
for nucleolin in chromatin dynamics by its histone chap- 
erone and co-remodelling activities (27). 

In this study, we analysed the association of nucleolin 
with rDNA and the consequences of nucleolin silencing on 
the transcription of pre-rRNA. ChlP-Seq analysis shows 
that nucleolin distribution along the rDNA coding and 
promoter regions is very similar to that of RPA116 and 
UBF. The silencing of nucleolin expression modifies the 
distribution of RNAPI, UBF, TTF-1 and TIP5 on the 
rDNA together with an increase of heterochromatin 



marks on the rDNA promoter and coding region. We 
propose a model in which nucleolin plays an important 
role in maintaining the euchromatin state of rDNA, 
allowing the efficient transcription elongation of rDNA 
by RNAPI. 

MATERIALS AND METHODS 

Cell lines and siRNA transfection 

HeLa cells were cultured in Dulbecco's modified Eagle 
medium (DMEM, Gibco) supplemented with 10% foetal 
bovine serum (Gibco) at 37°C in 5% C0 2 incubator. To 
generate the cell line stably expressing Flag-tagged 
nucleolin, HeLa cells were transfected with pSG5-Flag- 
nucleolin (encoding full-length human nucleolin with an 
amino-terminal Flag tag) and pcDNA™3.1 plasmids. 
Clones stably expressing Flag-tagged nucleolin identified 
by immunofluorescence were selected with hygromycin. 
HeLa cells (6 x 10 5 ) were transfected twice with siRNA 
(Eurogentec) using Lipofectamine 2000 (Invitrogen) as 
described previously (14). The siRNA sequences used are 
listed in Supplementary Table SI. Control siRNA was 
BLOCK-iT Alexa Fluor Red Fluorescent Oligo 
(Cat.No. 14750-100, Invitrogen). Proteins and RNAs 
were harvested at 96 h (for nucleolin siRNA) or 72 h 
after the first transfection (for TTF-1 siRNA) and 
analysed by western blot and reverse transcriptase-poly- 
merase chain reaction (RT-QPCR), respectively. 

Metabolic labelling and analysis of pre-rRNA 
transcription and processing 

Cells (4 xlO 5 ) were incubated at 37°C in 3 ml of DMEM 
without phosphate (Invitrogen) supplemented with 10% 
foetal calf serum, 25 mM HEPES (Invitrogen). After 2h, 
[ 32 P]Pi (PerkinElmer) was added to the medium 
(125 uCi/ml) for various times. At the end of the labelling, 
total RNA was extracted with RNeasy kit (Qiagen), 
separated by denaturing electrophoresis and the gel was 
dried before exposure. The quantity of RNA correspond- 
ing to an equal number of cells was loaded for each 
sample. Radioactivity was detected with a Fujifilm 
FLA-5100 scanner and the pictures were analysed with 
the Multigauge V3.0 software (Fujifilm). 

Transcription analysis by RT-QPCR 

Total RNA was prepared by TRIzol (Invitrogen) extrac- 
tion and digested with 20 U of RNase-free DNase I 
(Roche Diagnostics) for 1 h at 37°C. One hundred nano- 
grams of total RNA was reverse-transcribed using 
hexamer random primers and first-strand cDNA synthesis 
kit (Fermentas) and the synthesized cDNA was used for 
RT-QPCR using FastStart Universal SYBR Green 
Master (ROX) (Roche). Primer sequences used are listed 
in Supplementary Table S2. 

Immunofluorescence 

HeLa cells grown on glass coverslips were washed with 
phosphate buffered saline (PBS), fixed with 4% 
paraformaldehyde in PBS for lOmin at room temperature, 



Nucleic Acids Research, 2012, Vol 40, No. 19 9443 



incubated on ice for 20min and then permeabilized with 
0.1 % Triton X-100 in PBS for 2 x 5 min. After two washes 
in PBS, non-specific binding of antibodies was blocked 
with 10% foetal calf serum in PBS. After three washes 
with PBS, coverslips were incubated in primary antibodies 
at 37°C for 1 h. After two washes in PBS with 0.1% 
Triton X-100, coverslips were incubated with secondary 
antibodies coupled with Alexa dyes (A555 or A647). 
After two more washes in PBS with 0.1% Triton X-100, 
the coverslips were washed in PBS, rinsed in ddH 2 0 and 
briefly dipped in 100% ethanol. After a quick dry, 
coverslips were mounted with Fluoromount G containing 
200ng/ml DAPI. 

Chromatin immunoprecipitation and DNA 
methylation assays 

Chromatin immunoprecipitation (ChIP) assays were done 
as described earlier (28). HeLa cells were cross-linked for 
10 min with 1% formaldehyde in the culture medium at 
room temperature. Glycine (0.125 M) was added, and the 
cells were rocked for 5 min to stop the reaction. The cells 
were suspended in 200 ul sodium dodecyl sulphate (SDS) 
lysis buffer (1% SDS, lOmM EDTA, 50 mM Tris pH 8.1) 
and submitted to sonication to produce DNA fragments 
of 200-1000 bp in length. Chromatin was diluted 10-fold 
with ChIP dilution buffer (0.01% SDS, 1.1% Triton 
X-100, 1.2mM EDTA, 16.7mM Tris pH 8.1, 167mM 
NaCl), pre-cleared and immunoprecipitated with the re- 
spective antibodies. Precipitated DNA and protein 
complexes were reverse cross-linked and purified through 
phenol/chloroform extraction and ethanol precipitation. 
For RNase treatment, permeabilized HeLa cells were 
treated with 1 mg/ml RNase A for 5 min, cross-linked 
with formaldehyde and subjected to ChIP with 
anti-nucleolin antibody. 

For sequential ChIP (reChIP) assays, the specific DNA- 
protein complexes from the first immunoprecipitation 
were extracted by adding 25 JlxI of lOmM dithiothreitol, 
followed by incubation for 20 min at 37°C with vortexing 
every 5 min. The supernatants were pooled and diluted 10 
times with reChIP buffer (20 mM Tris-HCl pH 8.0, 2mM 
EDTA, 150mM NaCl, 0.1% Triton X-100), followed by 
the second immunoprecipitation. 

The purified DNAs were amplified by real-time PCR 
using StepOne Plus (Applied Biosystems) and FastStart 
Universal SYBR Green Master (ROX) (Roche). Primers 
for real-time PCR were designed or as reported (29) and 
are listed in Supplementary Table S2. 

For nucleolin ChlP-seq analysis, MOO ng of input DNA 
and DNA precipitated by nucleolin antibody were 
sequenced with the Illumina/Solexa 1G technique 
(GATC, Germany). The input DNA sequences 
(27 615 364 reads) were used for normalization of the 
nucleolin ChlP-seq data (22 381 018 reads). Unique reads 
from the nucleolin ChIP dataset (15 467484 reads) were 
mapped onto human rDNA reference sequence (NCBI 
accession number: HSU13369) using BWA (30,31). We 
optimized the quality threshold ( = 30) and mismatch 
penalty (— M) as 7 and -d as 12 to obtain the mapping 
of unique reads. The mapped nucleolin data and the 



control input data were then processed and filtered using 
tools and software algorithms from BioCOS Life Sciences, 
which uses data map quality-based filtering to eliminate 
repeated and poorly mapped reads. Furthermore, 
the mapped and filtered nucleolin data is normalized by 
the mapped and filtered input control data and plotted 
as shown on Figure 3A. ChlP-seq data of UBF and 
RPA116 (32) were downloaded from NCBI (NCBI sra 
SRR087746.sra (UBF HEK293T), SRR087747.sra 
(RPA116 HEK293T), SRR087753.sra (input HEK293T), 
SRR087754.sra (UBF K562) and SRR087755.sra (input 
K562) and analysed using the same methods (Supple- 
mentary Figure S3). 

To monitor CpG methylation of rDNA, the immuno- 
precipitated as well as the input DNA were digested with 
Hpall or Mspl, and the digested DNA along with an 
equal amount of the undigested DNA were amplified by 
real-time PCR using the H42.9 primers (Supplementary 
Table S2) specific for the human rDNA promoter region 
(29). The results are expressed as the ratio of methylated 
or unmethylated DNA to the respective input. 

Psoralen cross-linking and Southern blot 

HeLa cells grown in 60-mm Petri dishes (0.8 x 10 6 cells/ 
Petri) were placed in serum-free media (1.5 ml) just before 
cross-linking began. A 1/20 volume (75 (il) of 200(ig/ml 
trioxsalen (4,5 / ,8-trimethylpsoralen; Sigma) in methanol 
was added with gentle mixing, and after 5 min incubation, 
cells were irradiated on ice for 5 min with a 366-nm UV 
lamp (BlackRay model B-100A, 100W) placed at a 
distance of 6-7 cm. Cross-linking was repeated three 
more times, each time adding fresh Trioxsalen. The 
cross-linking procedure was performed in a dark room. 
Subsequently, the cells were washed with PBS and lysed 
with 0.5ml of lysis buffer (10 mM Tris, pH 7.5, 50 mM 
NaCl, 25 mM EDTA, 2% SDS). Proteinase K (1 mg/ml; 
Fermentas) was added, and the combined lysate was 
incubated overnight at 50° C, adding fresh proteinase K 
after the first 2 h. The genomic DNA was then de-proteini- 
zed twice with phenol/chloroform, ethanol precipitated 
and re-suspended in 250 pi TE pH 7.5, 0.1% SDS. 
Twenty micrograms of RNase A (Fermentas) was added 
and after 30 min at 37°C, 0.5 mg/ml Proteinase K was also 
added and incubation continued for another 1 h at 50°C. 
After two phenol/chloroform extractions, the genomic 
DNA was ethanol precipitated, re-suspended in 50 ul TE 
and quantified. About 10(ig of genomic DNA was 
digested overnight with BamHI and resolved on a 1% 
TAE agarose gel at 2Vcm _1 for 22 h in the absence of 
ethidium bromide (EtBr). The gel was subsequently 
EtBr-stained, photographed, its cross-links reversed by 
UV irradiation at 254 nm (4000 mJ- cm" 2 ) m a 
Stratalinker 1800 cross-linker (Agilent Technologies) and 
transferred onto a Biodyne B membrane (Pall). 
Hybridization was performed overnight at 65°C with 
labelled DNA probes in Southern Church Buffer (0.5 M 
sodium phosphate, pH 7.2, 1 mM EDTA pH 8.0, 
7% SDS, 1% BSA). Subsequently, membranes were 
washed with 6x saline sodium citrate (SSC), 2x SSC 
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and 0.1 x SSC and 0.1% SDS. Data were analysed using a 
FUJI Phosphorlmager and with MultiGauge software. 

Antibodies 

The following antibodies were used: a rabbit polyclonal 
antibody against human nucleolin (number 5567; de- 
veloped by our laboratory), anti-acetyl H4K12 (abl761; 
Abeam), anti-histone H3 (abl791; Abeam), anti-trimethyl 
H3K4 (ab8580; Abeam), anti-G9a (ab40542; Abeam), 
anti-TTF-1 (A302-361A; Bethyl), anti-TIP5 (CS-090-100; 
Diagenode), anti-dimethyl H3K9 (abl220; Abeam), 
anti-nucleolin monoclonal antibody (KAM-CP100; 
Stressgen), anti-RPA116 (a generous gift from Ingrid 
Grummt), anti-UBF (F-9) (sc-13125; Santa Cruz), 
anti-B23 antibody (abl0530; Abeam), anti-HDACl 
(ab46985; Abeam) and anti-Flag antibody (Sigma). 

RESULTS 

Depletion of nucleolin affects epigenetic marks of 
rDNA chromatin 

Previous experiments have shown that nucleolin is 
required for the efficient production of pre-rRNA by 
RNAPI, but the mechanisms implicated in this regulation 
have not been investigated (14,19,20). Transfection of 
HeLa cells with specific siRNAs against nucleolin 
efficiently reduced the levels of nucleolin protein 
(Figure 1A) and messenger RNA (Figure IB) and this 
was associated with a low level of pre-rRNA accumulation 
(Figure IB). Metabolic labelling with [ 32 P]Pi and northern 
blot analysis of pre-rRNA transcribed in nucleolin 
depleted cells confirm that nucleolin depletion severely 
impairs rRNA production without affecting the global 
maturation of pre-rRNA (Figure 1C and Supplementary 
Figure SI). By summing the intensities of all four major 
products, the loss of nucleolin causes a reduction of ~80% 
of rRNA accumulation compared to the control cells 
treated with a control siRNA. We then examined the mat- 
uration of the largest precursor (47S/45S) by comparing 
its intensity with the sum of all other maturation products 
and we observed that its maturation was only very mod- 
erately affected (decreased by ~10%) (Figure ID). The 
efficiency and rate of processing of the 32S species was 
not significantly affected (P > 0.05) in nucleolin-depleted 
cells (Figure IE and F) showing that 45S pre-rRNA pre- 
cursor is processed as in control cells. These data indicate 
that the low accumulation of pre-rRNA in nucleolin- 
depleted cells is not the consequence of an inhibition of 
rRNA processing that may affect ribosome biogenesis and 
transcription, but rather the result of a role of nucleolin on 
the transcription process itself. To better understand the 
mechanisms involved, we first investigated the conse- 
quences of nucleolin depletion on the epigenetic marks 
of rDNA chromatin (Figure 2). 

To determine whether this reduced level was associated 
with an increase in silencing by rDNA methylation, the 
level of methylation of a critical CpG dimer in the core 
promoter element of the rDNA was determined in control 
and siRNA- treated HeLa cells (Figure 2A). Nucleolin de- 
pletion did not affect the level of rDNA methylation at 



this site, suggesting that the decrease of pre-rRNA level is 
most likely not associated with increased rDNA methyla- 
tion. The level of active and inactive rRNA genes can be 
also determined by their susceptibility to the DNA 
cross-linking agent psoralen (33). To determine whether 
the decrease of pre-rRNA synthesis was associated with 
the alteration of the chromatin state of rDNA, we 
performed this psoralen cross-linking assay. We found 
that nucleolin depletion has no significant effect on the 
proportion of genes that are efficiently cross-linked by 
psoralen, suggesting that the ratio of active to inactive 
rRNA genes is not changed in nucleolin-depleted cells 
(Supplementary Figure S2). Then, we examined the 
histone post-translational modifications on rDNA 
promoter and coding regions (Figure 2B). In agreement 
with the lower accumulation of pre-rRNA molecules in 
nucleolin-depleted cells, we observed an important 
decrease in H4K12Ac and H3K4me3, two marks usually 
associated with transcriptionally active rDNA repeats, 
whereas the level of H3K9me2 increased significantly. 

As nucleolin is required for rRNA transcription, we 
could also expect that over-expression of nucleolin 
affects the chromatin state of rDNA and pre-rRNA tran- 
scription. Indeed, in HeLa cells that stably over-express 
Flag-tagged nucleolin (Figure 2C), we observed a marked 
increase of pre-rRNA level (Figure 2D). This was 
associated with a strong increase in the activation marks 
H3K4me3 and H4K12Ac and a decrease in the silencing 
mark H3K9me2 in rRNA genes (Figure 2E). Altogether, 
these ChIP data indicate that the depletion and over- 
expression of nucleolin are associated with opposing alter- 
ations of histone modification states of rDNA but not 
with a change of the proportion of active genes. 

Nucleolin is associated with rDNA 

These alterations of post-translational modifications of 
rDNA chromatin upon nucleolin silencing or over- 
expression could be indirect or the result of nucleolin 
interaction with rDNA chromatin. To characterize these 
interactions, we performed ChlP-Seq experiments 
with nucleolin antibody. ChIP experiments from two in- 
dependent experiments were analysed. The total number 
of reads from nucleolin ChlP-seq (22 381 018 reads) and 
from input DNA (27 615 364 reads) were mapped on the 
genome, and only unique reads (15 467484) were mapped 
on rDNA (NCBI accession number: HSU13369). The 
mapped nucleolin data and the input starting material 
were then processed and filtered as described in the 
'Materials and Methods' section, and the nucleolin 
dataset was then normalized by the starting input 
material and plotted (Figure 3 A). These results show 
that nucleolin is particularly enriched in the coding 
region of rDNA, moderately enriched in the promoter 
region and insignificantly present in the IGS. Sequence 
reads that map to multiple locations were masked 
during bioinformatic analysis, which can leave gaps in 
regions of rDNA that are rich in repetitive sequences; 
this is particularly observed in the 18S region (star in 
Figure 3A). The distribution of nucleolin on rDNA was 
compared to the distribution of UBF and RNAPI 
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Figure 1. Effect of nucleolin on rRNA transcription and processing. (A) Inhibition of nucleolin expression by siRNA. HeLa cells were transfected 
with control (BLOCK-iT™) or nucleolin siRNA. Proteins (nucleolin and B23) were analysed by western blot 4 days after transfection. **P<0.01, 
Student's /-test was done between HeLa cells transfected with control siRNA and nucleolin siRNA. (B) Levels of nucleolin messenger RNA and 
pre-rRNA determined by RT-QPCR as described in the 'Materials and Methods' section. The data represent the average of three different inde- 
pendent experiments. **P<0.01, Student's /-test was done between HeLa cells transfected with control siRNA and nucleolin siRNA. (C) Metabolic 
labelling of total RNA from HeLa cells transfected with control siRNA or nucleolin siRNA. HeLa cells were incubated in phosphate-free DMEM 
for 2h and [ 32 P]orthophosphate was added to the medium (125fiCi/ml) for the indicated times. The 32 P-labelled RNAs corresponding to an equal 
number of cells were resolved on 1% MOPS-formaldehyde gels and exposed on a Fujifilm FLA-5100 scanner. (D) Quantification of the relative 45S 
processing in HeLa cells transfected with control and nucleolin siRNA as in (A). For each pulse time, the 45S processing was calculated as the sum 
of the intensities of processed products 32S, 28S and 18S, divided by the intensity of the 45S precursor. Data were normalized to the control cells. 
Data were from two independent experiments. (E) Quantification of the relative 32S processing in HeLa cells transfected with control and nucleolin 
siRNA as in (A). For each pulse time, the 32S processing was calculated as the intensity of processed product 28S divided by the intensity of the 32S 
precursor. Data were normalized to the control cells. Data were from two independent experiments. (F) Quantification of the rate of 32S rRNA 
processing in HeLa cells transfected with control and nucleolin siRNA. The 32S rRNAs were normalized to the 45S rRNAs. Data were from two 
independent experiments. Student's /-test was done between HeLa cells transfected with control siRNA and nucleolin siRNA, and no significant 
difference was detected (panels D, E and F). 
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Figure 2. Effect of nucleolin on epigenetic marks. (A) Depletion of nucleolin does not affect rDNA promoter methylation. DNA was isolated from 
control or siRNA-treated HeLa cells and digested with Mspl (methylation-insensitive) or Hpall (methylation-sensitive). The human rDNA promoter 
was then amplified with primers flanking the CCGG site at —9. If the CpG dimer was methylated (resistant to cleavage with Hpall), amplification of 
rDNA from —57 to +43 yield a fragment of 100 nt, whereas unmethylated DNA would be cleaved and yield no PCR product. Student's /-test was 
done between HeLa cells transfected with control and nucleolin siRNA, no significant difference was detected. (B) Histone post-translational 
modification marks of rDNA are affected after the depletion of nucleolin. QPCR amplification of either the S'-end of the rRNA gene (H42.9) or 
the coding region (H8) after ChIP shows the level of H3, H3K4me3, H3K9me2 and H4K12Ac on different regions of rDNA. The ratio of rDNA 
occupancy between nucleolin and control siRNA transfected cells was calculated for each histone mark. (C) Over-expression of Flag-nucleolin. Total 
cell extract from the HeLa cell line stably expressing Flag-nucleolin was used for western blot analysis with an anti-nucleolin and anti-tubulin 
antibody, and the quantification was done using Image J. (D) The level of pre-rRNA in cells over-expressing nucleolin was determined by RT-QPCR 
with primers specific for the 5 ; -ETS. (E) Over-expression of nucleolin affects histone post-translational modifications on rDNA. ChIP experiments 
shows the level of H3K4me3, H3K9me2 and H4K12Ac on 5'-end of the gene (H42.9) and coding regions (H8) of rDNA in cells that stably 
over-express Flag-nucleolin. For all experiments, *0.01 < P < 0.05, **P<0.01, Student's /-test was done between control HeLa cells and HeLa 
cells transfected with nucleolin siRNA, or cells over-expressing Flag-tagged nucleolin. 
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Figure 3. Nucleolin is associated with rDNA. (A) ChlP-seq mapping of nucleolin binding sites throughout the rDNA locus. Data were analysed as 
explained in the 'Materials and Methods' section, and the nucleolin tag density normalized to the starting material was plotted. As reads that map to 
multiple locations were masked during the bioinformatic analysis, regions rich in repetitive sequences can leave gaps as shown in the 18S region 
(star). (B) Schematic representation of a human rDNA repeat. The positions of QPCR amplicons in ChIP assays are indicated with solid bars. 
(C) QPCR analysis of the enrichment of nucleolin on rDNA. Immuno-precipitated DNA was analysed by QPCR using sets of primers indicated in 
(B). The percentage of DNA immuno-precipitated with anti-nucleolin antibody was calculated relative to the ChIP input DNA. Values are 
means ± SD (standard deviation) derived from three independent ChIP experiments, each tested by at least three independent QPCRs. (D) HeLa 
cells stably transfected with control Flag-EV (white bars) or with Flag-tagged nucleolin (grey bars) were used for ChIP experiments using anti-flag 
antibody. The precipitated DNA was analysed by QPCR as previously. The percentage of DNA immuno-precipitated with anti-Flag antibody was 
calculated relative to the ChIP input DNA as in Panel C. 
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(RPA116) that have recently been published (32) 
(Supplementary Figure S3). Interestingly, the global 
pattern of nucleolin distribution on rDNA is very 
similar to the distribution of UBF and Poll. To validate 
these ChlP-seq data, we performed several QPCR 
analyses using primers that span the human rDNA 
repeat (29) (Figure 3B and C). As expected from the 
ChlP-Seq data, nucleolin was preferentially found 
associated with the rDNA-coding regions, including the 
18S rRNA, which does not show any unique reads in 
the ChlP-seq analysis (star in Figure 3A). The interaction 
of nucleolin with rDNA chromatin was independent of 
RNA since RNase treatment before the ChIP procedure 
does not abrogate the nucleolin with the rDNA locus 
(Supplementary Figure S4). Furthermore, when the 
ChIP experiment was performed with a cell line that 
stably expresses Flag-nucleolin (Figure 3D), we observed 
the similar enrichment of nucleolin in the rDNA-coding 
sequence. 

To determine if the decreased accumulation of 
pre-rRNA in nucleolin-depleted cells was linked to a modi- 
fication of the interaction between RNAPI and UBF with 
the rRNA genes, we performed ChIP experiments using 
the antibodies to the RNAPI subunit (RPA116) and for 
UBF in nucleolin-depleted cells. Nucleolin depletion leads 
to a decrease in UBF on the promoter and coding regions 
of rDNA (Figure 4 A and C), which is in agreement with 
the decrease in rRNA transcription in nucleolin-depleted 
cells. Interestingly, nucleolin depletion leads also to an 
increase of Poll loading at the beginning of coding 
region of rDNA but to a significant decrease at the 
3 / -end of the gene (Figure 4B and C). This pattern of 
RNAPI loading after nucleolin depletion is reminiscent 
of what happens after inhibition of transcription elong- 
ation by actinomycin D (34), suggesting that nucleolin 
could be required for efficient RNAPI elongation. 

Nucleolin is associated with active rDNA repeats 

We then determined the epigenetic state of the rDNA 
chromatin associated with nucleolin (Figure 5). Active 
and silent rDNA can be distinguished by DNA methyla- 
tion. Active rDNA is hypo-methylated, while silent rDNA 
is hyper-methylated (6). In human cells, methylation at a 
single site in the promoter of the rRNA gene correlates 
with the repression of the promoter activity (35). We per- 
formed a ChlP-CHOP experiment to look at rDNA 
methylation status at the —9 CpG site of the promoter 
sequences bound to nucleolin (Figure 5A). As a control, 
we did the same experiment with UBF, which is almost 
exclusively associated with unmethylated rDNA pro- 
moters (36). Interestingly, nucleolin was found to be also 
preferentially associated with unmethylated rDNA pro- 
moters, which correspond to the transcriptionally active 
rRNA genes. Re-ChIP experiments indicate that the het- 
erochromatin mark H3K9me2 in the nucleolin-associated 
chromatin seems absent within the rDNA coding region 
(Figure 5B) but is highly enriched in the promoter 
sequence (Figure 5C). In contrast, H4K12Ac and 
H3K4me3 euchromatin marks are indeed associated 
with nucleolin in the rDNA coding sequence (Figure 5B) 



and at the 5'-end of the rRNA gene (H42.9 primers) 
(Figure 5C), which is in agreement with an association 
of nucleolin with the transcriptionally active copies of 
rDNA. 

Nucleolin affects TTF-1 interaction with TO 

TTF-1, which binds the promoter-proximal terminator 
TO, plays a key role in the recruitment of the nucleolar 
remodelling complex NoRC and the silencing of rDNA 
(5), although earlier reports indicated that TTF-1 was 
required for the transcription by RNAPI of chromatinized 
DNA templates (37). In order to determine whether 
nucleolin plays any role in this process, we studied 
whether the interaction of nucleolin with ribosomal 
chromatin had any influence on the binding of TTF- 1 to 
the TO site (Figure 6). We found that nucleolin depletion 
leads to a two-fold increase in TTF-1 interaction with TO 
(Figure 6B). As the binding of TTF-1 is required for the 
recruitment of NoRC and HDACs, we also performed 
ChIP experiments with TIP5 and HDAC1 after nucleolin 
depletion (Figure 6C). We also observed a slight increase 
in the recruitment of TIP5 and HDAC1, while the inter- 
action of G9a was not changed. In a reverse experiment, 
we observed that the inhibition of expression of TTF1 
leads to increased association of nucleolin with TO 
(2-fold) and H42.9 (4-fold) (Figure 6D). Altogether, 
these data show that the presence of nucleolin on rDNA 
chromatin inhibits the binding of TTF-1 within the TO 
region and vice versa and thus protects rRNA genes 
from TTF-1 -mediated silencing of transcription. 

DISCUSSION 

Over the past 20 years, many in vitro and in vivo studies 
have implicated nucleolin, one of the major nucleolar 
proteins, in the production of rRNAs by RNAPI tran- 
scription (14,15,19-21,38) without providing many mech- 
anistic details on how nucleolin could participate in the 
production of rRNA. Previous works have shown that in 
HeLa cells the accumulation of 45 S could be affected by 
the rate of pre-rRNA processing (39,40). Since nucleolin 
interacts specifically with pre-rRNA (41-47) and has been 
involved in the first processing step of pre-RNA in vitro 
(16), it was tempting to explain the low accumulation of 
pre-rRNA in nucleolin depleted cells by an indirect effect 
of nucleolin on pre-rRNA processing. However, by meta- 
bolic labelling or northern blot we could not detect major 
changes in the processing pathways of pre-rRNA or in the 
efficiency of this processing (Figure 1 and Supplementary 
Figure SI) that could explain the strong reduction of 45S 
accumulation. These data are also in agreement with our 
previous analysis of nucleolin knockout in chicken DT40 
cells (14,20). One possible explanation is that the low level 
of nucleolin that remains in nucleolin-depleted HeLa cells 
is sufficient to support normal pre-rRNA processing, 
while it is affecting very strongly pre-rRNA accumulation 
through its transcription. Indeed, the accumulation of 
pre-rRNA is very sensitive to the level of expression of 
nucleolin (20). We have seen the same effect of nucleolin 
depletion on the level of pre-rRNA not only in HeLa cells 
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but also in human primary fibroblast (14) and in chicken 
DT40 cells (20), showing that what we describe in this 
article is not a particularity of HeLa cells. 

Recent reports indicate that nucleolin is a histone chap- 
erone with a FACT-like (FAcilitates Chromatin Tran- 
scription) activity that facilitates chromatin remodelling 
and histone dynamics (26,27), suggesting that nucleolin 
could regulate rRNA transcription at the level of chroma- 
tin accessibility and transcription elongation. 

In this report, we describe experiments that support the 
role of nucleolin in the maintenance of the euchromatin 
state of rRNA, which allows efficient RNAPI transcrip- 
tion elongation. We show that nucleolin is preferentially 
associated with unmethylated rDNA and with epigenetic 
marks characteristic of transcriptionally active genes, 
which is in agreement with earlier reports indicating that 
nucleolin is required for rRNA transcription. In 
Arabidopsis thaliana, disruption of the AtNUC-Ll gene 
induces the loss of DNA methylation without affecting 
histone epigenetic marks at rRNA genes (48). However, 
we found that, in human cells, the depletion of nucleolin 
strongly inhibits pre-rRNA transcription and affects the 
post-translational modifications of histones associated 
with the rDNA promoter (Figure 2B and E) but does 
not modify the level of CpG methylation (Figure 2A) in 
the promoter region. This suggests that in mammalian 
cells, increased CpG methylation is not required for the 
inhibition of rRNA transcription in response to nucleolin 
depletion. Of note, methylation-independent silencing of 
rRNA genes was also observed upon depletion of UBF 
(36). Interestingly, UBF depletion leads to an increase in 
the number of rRNA genes in an inactive condensed state 
without reducing the net rRNA synthesis as transcription 
from remaining active genes is increased (36). Upon 
nucleolin depletion, we observed a substantial decrease 
in UBF loading on rRNA gene (Figure 4); in contrast to 
Sanij et al. (2009) report, this is correlated with a strong 
decrease of rRNA accumulation, as nucleolin is probably 
required for the transcription activation of the genes that 
remained actives in UBF-depleted cells (36). 

Nucleolin depletion does not affect the proportion of 
active rRNA genes detected by the psoralen cross-linking 
experiment (Supplementary Figure S2). It should be noted 
that there is not always a direct correlation between the 
level of psoralen cross-linking and the presence of hetero- 
chromatin marks or the level of RNAPI transcription. 
For instance, when DNA methylation is lost, 50% of 
rRNA genes still show an inactive state by psoralen (49). 
In addition, UBF depletion leads to a decrease in the 
active ribosomal genes in psoralen cross-linking experi- 
ments (36); however, rRNA synthesis does not decrease 
as the transcription from the active genes is increased. This 
result may suggest that, in nucleolin-depleted cells, the 
lower level of pre-rRNA accumulation is rather the con- 
sequence of a lower rate of transcription of the active 
genes rather than an increase in the number of inactive 
genes. Indeed, our data showed that nucleolin depletion 
leads to an increase in RNAPI at the 5'-end of the coding 
region of the rDNA (Figure 4B and C), as previously 
described for the effect of actinomycin D on transcription 
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elongation (34), which is consistent with a role for 
nucleolin in RNAPI elongation. 

In our ChlP-seq analysis, we found a strong enrichment 
of nucleolin at the rDNA coding region with little if any 
enrichment in the IGS. Interestingly, this distribution 
looks very similar to that of UBF, another factor 
required for rRNA transcription (32) and to that of 
RNAPI subunit RPA116. The association of nucleolin 
with rDNA does not result from its association with 
pre-rRNA precursor, as RNase treatment does not alter 
the distribution of nucleolin on rDNA (Supplementary 
Figure S4). The presence of nucleolin on the coding 
region is associated with the euchromatin marks 
H4K12Ac and H3K4me3 (Figure 5B and C). Although 
H3K9me2 is usually associated with heterochromatin for- 
mation, there are now several reports that demonstrate 
that components of heterochromatin such as H3K9me2 
play additional role in establishing chromatin structure 



of actively transcribed genes. This is true for PolII genes 
(50) and also for Pol I genes (51) indicating that the asso- 
ciation of nucleolin with H3K9me2 at the 5'-end of the 
gene may be another sign of the link between active genes 
and nucleolin. Upon nucleolin depletion, the level of 
H4K12Ac and H3K4me3 on rDNA decreased while the 
level of H3K9me2 increased. On the other hand, nucleolin 
over-expression caused an increase in H3K4me3 and 
H4K12Ac and a decrease in H3K9me2, suggesting that 
the presence of nucleolin is required for the modulation 
of histone modifications of rDNA associated with rRNA 
transcription. 

The binding of TTF-1 to TO is an important step for the 
silencing of rRNA as this allows the recruitment of the 
nucleolar remodelling complex NoRC (3,5,6), which then 
recruits HDAC and HMT to establish a heterochromatin 
structure of rRNA (1-3). Interestingly, we found that after 
nucleolin depletion, the interaction of TTF-1 with TO 
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Figure 7. Model depicting the role of nucleolin in rRNA transcription. (A) Nucleolin is enriched in the promoter and coding regions of rDNA and 
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TIP5, then recruits HDAC and HMT to the rDNA promoter to establish heterochromatin structure thus inhibit rRNA transcription. 



increases, with a concomitant increase of TIP5 on rDNA 
promoter (Figure 6B and C). This is also accompanied by 
recruitment of HDAC 1 (Figure 6C), in agreement with the 
lower level of H4K12Ac that we observed (Figure 2B). 
NoRC function requires the interaction of its TIP5 
subunit with 1 50-250 nt pRNAs transcribed from a 
spacer promoter located ~2kb upstream of the major 
45S pre-rRNA promoter (52). These pRNAs also seem 
to interact with TO leading to displacement of TTF-1 
and to the recruitment of DNA methyltransferase 
DNMT3b to the promoter (13). After nucleolin depletion, 
even if we observed higher recruitment of TIP 5, we could 
not detect any difference in the —9 bp CpG methylation, 
suggesting that this mechanism for heterochromatin for- 
mation and rDNA silencing mediated by NoRC does not 
entirely take place in nucleolin-depleted cells. 

In addition, nucleolin-depleted cells show a striking 
reorganization of nucleoli structure (14,20,53), and 
several reports show that nucleolin is implicated in chro- 
matin condensation (54), chromatin loop organization 
(55,56) and histone dynamics in vivo (27). Therefore, 



nucleolin could also be implicated in the regulation of 
chromatin accessibility at a relatively large scale in 
addition to its effect on single nucleosome co-remodelling 
and destabilizing activity (26,57). 

Altogether, this study provides new information to 
explain the role of nucleolin in regulation of RNAPI tran- 
scription (Figure 7). Nucleolin depletion results in an 
increase of RNAPI at the 5'-end of the coding region of 
the rDNA and a decrease of UBF along the rDNA. The 
association of nucleolin with the promoter region inter- 
feres with the binding of TTF-1 to TO, thereby inhibiting 
the recruitment of NoRC. Nucleolin depletion causes an 
increase of TTF-1 on TO, allowing a better recruitment of 
NoRC and then the establishment of a heterochromatin 
state contributing to the inhibition of RNAPI 
transcription. 



SUPPLEMENTARY DATA 
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